Comment on Perez et al PRX 2, 041005 (2012) 



A. Beresnyak 
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Recently Perez et al [lt| wrote on the spectral slope of MHD turbulence claiming that it is consistent 
with -3/2. This work contains a number of errors, factual inaccuracies and puzzling methods in 
the interpretation of numerics. We argue that the numerical evidence that the authors presented 
is actually against -3/2 slope and the assumption of universal alignment, opposite to what they 
claim. Perhaps the most puzzling is the Fig. 8 that claims that the authors measured the inertial 
range length and that it is consistent with Re^^^ scaling. At a close inspection it appears that the 
datapoints are not the result of a measurement, but rather were calculated by a formula. 
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Currently there were conflicting reports regarding the 
spectral slope of MHD turbulence. While Beresnyak re- 
ports numerical data as roughly consistent with —5/3 
slope and inconsistent with —3/2 0,01, Perez et al [l| 
group claims it is consistent with —3/2. Both groups 
do driven simulations of strong turbulence in the strong 
mean field limit by using pseudospectral code. For dis- 
sipation scheme Beresnyak used hyperdiffusion as well 
as normal diffusion, while Perez et al used only normal 
diffusion. 

Recently Perez et al published a detailed paper on their 
latest high-resolution simulations with long evolution in 
time, and also claimed that Beresnyak's result is dom- 
inated by numerical error. The paper is the subject of 
this comment. 

In this comment we discuss several major issues with 
Perez et al interpretation of numerics, including: a) the 
issue of Fig. 8, b) the claim of the universality of align- 
ment, c) the claim of —3/2 slope, d) the claim of the 
"fake" numerical convergence. 

The issue of Fig. 8 — On Fig. 8 the Authors claimed 
that they measured the length of the inertial range and 
that its scaling with Reynolds number is consistent with 
Boldyrev's model [4, 5]. On close inspection of Fig. 8 it 
is evident that the Authors arbitrarily took a constant 
A:^ = 4 as the beginning of the inertial range and the 
constant k±ri = 0.1 as the end of it. The length of the in- 
ertial range that the Authors "measured" and presented 
in the lower panel of Fig. 8 was, therefore, calculated by 
a formula lo/ld — 4077^^, where 77 is defined by Eq. 8, 
so la /Id = 40e2/9A-i/9j/-2/3^ All quantities that enter 
this formula are the parameters of the simulations, so 
this "measurement " does not constitute a real measure- 
ment. Also with e^/^ vary only by a factor of 0.99 and 
A^/^ vary by a factor of 0.99, the main dependence of 
lo/ld is, automatically, v""^/^ ~ i?e^/^, which was con- 
firmed on the lower panel and incorrectly attributed to 
the correspondence with Boldyrev's model. If one does a 
real measurement assuming, e.g., that the inertial range 
length is the greatest stretch of the spectra, deviating 
from —3/2 law by no more than 5%, then the scaling 
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FIG. 1. The inertial range length reported by Perez et 
al (crosses), which they calculated by a formula and the 
overplotted results of the real measurement of such length 
(squares), assuming that the inertial range length is the great- 
est stretch of the spectra, deviating from —3/2 law by no more 
than 5%. 



of such inertial range length is grossly inconsistent with 
Re^/^, see Fig 1 of this comment. The same is true if 
some other criterion is used, e.g., using deviation of 10%. 

We conclude that the Author's argumentation that nu- 
merics exhibit Re^l^ scaling of the inertial range length 
is errorneous and void. Instead, the measurements of the 
inertial range length show a non-power law dependence 
on Re. 

As a side-note, the length of the —3/2 range that Perez 
et al simulations seem to saturate to with increasing Re, 
8 10, is consistent with the length of the —3/2 range 
published in 0, 0] ■ 

The claim of the -3/2 slope — Perez et al simula- 
tions feature rather long evolution in time, 100 dynamical 
times on the outer scale. These simulations, therefore, 
will produce measurements with significantly lower sta- 
tistical error than that of others, which typically used 
10-20 outer dynamical timef0. For the purpose of this 
comment wc will only assume that the statistical error in 



^ It have to be noted that smaller scales have much better statistics 
because they have smaller timescales and much more realizations 
in the simulation box, which is why 10-20 outer dynamical times 
is usually considered adequate for the turbulence scaling studies. 
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Perez et al is lower than that in previous hydrodynamic or 
MHD simulations, such as The scaling conver- 

gence argument which has been used in the above studies 
require that all curves from different simulations collapse 
into a single curve starting with some scale / and on all 
smaller scales, including dissipation scales. If there is a 
convergence in the inertial range but there is no conver- 
gence on the viscous scales, this would mean that there 
is some source of systematic error on grid scale that ac- 
tually depends on the resolution of the experiment, such 
as the numerical equations on grid scale are actually dif- 
ferent for 1024^ and 2048^. The author's Fig 8 does show 
such a lack of convergence within k±ri = 0.08 1.0. The 
authors, however, neither mention their non-convergence 
on these scales, nor they claim any such strange system- 
atic error. This error has never been reported before and 
[1, H, IB-Ill always report that the curves collapse onto the 
same curve starting with some scale I and on all smaller 
scales. To the best of our knowledge such systematic er- 
ror simply does not exist and true convergence should 
always be in the above conventional sense. 

As a curious excersize one can also check that if the 
perpendicular spectrum converges on the inertial range 
but does not converge on viscous scales, the other type 
of spectra, e.g. one-dimensional spectrum will not con- 
verge at alH. So, in this case the fact of the convergence 
hangs on the type of spectrum one uses, which is highly 
doubtful. In other words. Fig 8. does not demonstrate 
convergence of spectra with -3/2 scaling, on the contrary, 
it demonstrates the lack of such convergence. 

Also, it is quite puzzling that while in sections 2-3 the 
authors spent a great deal of time explaining why simu- 
lations with reduced parallel resolution are bad, they use 
the very same reduced parallel resolution simulations for 
convergence study, even though they have data from full 
parallel resolution simulations? Another puzzling fea- 
ture of Fig. 8 is the high level of fluctuations of RB3 
spectrum. All other spectra, presented in this paper are 
very smooth. The authors could have used a smaller time 
evolution for this spectrum but they never say so in the 
paper. 

The claim of the universality of alignment — Similar 
convergence arguments can be used to claim a universal 
scaling of some other quantity. The authors were trying 



^ One-dimensional spectrum is defined as a power spectrum of a 
quantity sampled along a certain line and then averaged over all 
lines and directions. This type of spectrum is of interest because 
it is produced by the solar wind measurements from satellites. 
The mathematical relation between isotropic or perpendicular 
spectrum E{k) and the one dimensional spectrum is P{k) = 
J'J"i?(fcl)^^ [^. From this expression it is clear that if two 
E\{k) and E2{k) are coincident in the inertial range but different 
on viscous scale, the corresponding Pi(fc) and P2(fe) are different 
everywhere. Note that llOll confusingly compare simulated E(k) 
and the measured P(k). 



to support universality of so-called dynamic alignment 
9{l). Again their figure Fig. 7 does not show convergence 
and having very low level of statistical error and no source 
of systematic error the authors should have been expect- 
ing rescaled dynamic alignment curves to collapse onto 
a single curve. They, however, do not, quite consistently 
with [3] . Contrary to the numerical evidence the authors 
still claim that 9{l) scalings are universal. 

The logical loop in the alignment universality claim - 
The Authors bring several qualitative arguments in favor 
of universality of their alignment measure 9{l), however, 
before they do a scaling study to demonstrate such uni- 
versality, they say that they rather postulate it and re- 
define outer scale A so that 9 {I) will fit better between 
simulations with different resolutions and after that say 
that alignment is "remarkably stable" . This argumenta- 
tion is, obviosly, a logical loop. But even after the Au- 
thors fitted their A, the 9{l){h/rif/'^ stiU didn't converge 
as noted above (Fig. 7). 

Defining outer scale by alignment is certainly un- 
usual. Conventionally defined outer scale, L = 
{3tt/AE) k~^E{k)dk is a robust measure that shows 
only minuscule dependence on resolution, as long as the 
driving procedure is unchanged. The Author's A, defined 
by alignment, change considerably, however. The Au- 
thors did not comment of whether they explicitly changed 
driving procedure between three simulations RBla-3a or 
why A could be different in each of these simulations, 
especially considering that e is virtually constant. 

Obviously, the A, defined by the Authors, depends on 
resolution, because alignment is not universal. 

The claim of the "fake" numerical convergence — In 
their Appendix the Authors claim that the convergence 
on viscous scales observed in [1, [3] is a "fake" or a nu- 
merical artefact convergence which is due to numerical 
error. This is fairly puzzling claim, considering the scal- 
ing study is a traditional method which has been used 
for years, e.g. Also it is quite puzzling because the 

convergence on viscous scales from different simulations 
when one compensates spectrum by Kolmogorov scaling 
and plots it against kr]4i simply means that the velocity 
perturbation on viscous scale is proportional to the Kol- 
mogorov velocitjH Vrj4i = (ve)^/'^. Why would numerical 
error conspire in such a way as to produce such specific 
dependence of perturbation amplitude on e and vl Why 
such a puzzling phenomenon has never been reported be- 
fore? 

In fact the scaling study is relatively unaffected by 
what happens on small scales, because it is just a rescal- 



Indeed, if the spectrum E{k) is compensated by e 2/3^5/3 g^^j 
converge to some constant C at krm equal to, say, 0.2, then 
the velocity perturbation 6v = ^ E{k)k on this scale will be 
determined by C^/'^0.2-^/^{uey/^. 
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ing argument [s*! . The difference between rescaled quanti- 
ties on the viscous scale have very httle to do with what 
happens on viscous scales, but has to do with the dif- 
ference in resolutions and the suggested rescaling slope. 
For example, taking 1024^ and 2048'^ simulations the lat- 
ter spectrum has to be shifted left by a factor of 2 and 
the spectrum has to be shifted up by a factor of 2^/^. 
Why would numerical error know that it has to produce 
a factor of 2^/^ less numerical noise on the grid scale in 
2048^ compared to 1024^? Why would numerical error 
on the grid scale be so different in these two simulations 
and why would it be different by this precise factor? 

The Authors also claim that choosing Nrjj^i — const 
resolution criteria, i.e. the conventional criteria that peo- 
ple used before, will result in some "numerical conver- 
gence" , which "should not be confused with the conver- 
gence to the physical solution". It seems that the Au- 
thors are unaware of the standard resolution studies that 
were performed to confirm the convergence of spectra by 
increasing resolution and keeping all other parameters 
constant. In a proper scaling study each individual spec- 
trum from each simulation is converged to the physical 
solution within some error. Therefore, convergence of nu- 
merical spectra means absolute convergence. Ironically, 
the Authors themselves use Nrj/^i = const resolution cri- 
teria and never comment on why they did not believe 
their own model, which they claim their data supported, 
and use Nrjn = const instead of Nrj = const. 

The Authors also claimed that their simulations are 
better resolved than 0, Q • In doing so they referred to 
the simulations R8 and R9 in [sj which are the only simu- 
lations presented in that paper, which are under-resolved 
by conventional criteria. They do not mention other sim- 
ulation groups, such as R6-7, which are quite similar in 
numerical setup to what the Authors presented in their 
resolution study in Fig. 8, i.e. RBla-RB3a. Indeed R6 is 
384 X 10242, while RB2a is 256 x 1024^, R7 is 768 • 2048^, 
while RB3 is 512 • 2048^. Moreover, R6-7 have lower 
Reynolds numbers and R6-7 has fcniax'74i — 0.95, while 
in RBla-3a kmaxVn = R6-7 are better resolved 

by a factor of 1.19, and, at the same time, have higher 
parallel resolution than RBla-3a. In other words, when 
comparing what the authors presented in their resolu- 
tion study on Fig. 8 and the resolution study of R6-7 
in the latter are better resolved in both parallel and 
perpendicular directions. The Authors, however, never 
mention this and made it sound like the simulations in 
Q are under- resolved compared to theirs. 

The Authors also make vague analogies between sim- 
ulated truncated Euler equations that shows "thermal- 
ized" tail and the driven hyperviscous simulations. At 
a close inspection these analogies do not hold. Indeed, 
the tail is essentially non-stationary, as the energy is 
being dumped from the turbulent cascade into this tail. 
One can only speak of a certain snapshot of such sim- 
ulation whose spectra might be similar to the spectrum 




FIG. 2. The spectra for imbalanced simulation reported in 
[l^ . compensated by the factors corresponding to [l3| model 
and [11 model. 



of the very high order hyperviscosity with gigantic bot- 
teneck bump. Out of simulations presented in 2, 3] only 
one has a visible bottleneck bump and its relative ampli- 
tude, 0.22 is even lower than the amplitude of bottleneck 
in hydrodynamics with normal viscosity, 0.31. Then, 
why hydrodynamic viscous scales are not similar to the 
"thermalized" tail and our viscous scales are? Further- 
more, the Authors speculate that the short —5/3 range 
observed [2] is the result of an anti-bottleneck effect. But 
that would require that the anti-bottleneck effect has an 
amplitude comparable to the amplitude of the bottle- 
neck effect itself. This has never been observed and the 
anti-bottleneck effect that the authors call "pseudodissi- 
pation" is always much weaker than bottleneck effect, is 
not noticeable in viscous hydro with its 0.31 amplitude 
and is not supposed to be visible in our simulation R3 
with its 0.22 amplitude of the bump. Needless to say, 
no scaling study argument is applicable to the "ther- 
malized" tail, which is not even statistically stationary 
and this tail is not actually thermalized [ll'|. As to the 
numerically resolved statistically stationary simulations 
with hyperviscosity, the scaling study argument is still 
well-applicable [3]. 

Furthermore, the authors claim that the alignment is 
"partially lost" in the "thermalized" region. This is, 
again, puzzling, because the alignment is not lost, i.e. 
its value does not equal to the value corresponding to 
the random vectors of v, B, but it just flattens out, i.e. 
becomes independent on scale as was clearly presented in 

a- 

Finally, on their last plot of the Appendix the Authors 
make a convergence study between simulations with dif- 
ferent geometries of the grid cell, RB2c and RB3a. They 
don't obtain any convergence with either -3/2 or -5/3. 
What can be derived from such a study is unclear, be- 
cause the scaling study argument 0, 0, 0] simply does 
not work in this case and one is not supposed to do con- 
vergence study between such simulations. The Authors, 
nevertheless, claim that this lack of convergence somehow 
supports their argumentation. 

Imbalance spectra - The Authors refer to previous stud- 
ies of imbalanced case as producing "conflicting results" 
and claim that their data further support the -3/2 spec- 
tral slope. However, as we see previously, even the bal- 
anced slope is not agreed upon and neither of the Au- 
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thor's spectra shows clear -3/2 slope. In this situation 
a rigorous resolution study is, again, necessary to con- 
firm or reject models. In particular, one wants to com- 
pensate spectrum by £;(fc)fc±(e±)-2/3L-i/6fc3/2 ^ jg 



correct or by E{k)^{e^)~^/^{e'^)'^/^k^/^ if [iJI is correct. 
Fig. 2 of this comment shows such comparison from a low- 
imbalance simulation previously reported in 1^. As we 
see the model is grossly inconsistent with numerics. 
Unfortunately, the Authors neither mention this nor they 
discuss how the predictions of should be modified in 
order to be consistent with numerics. Furthermore, citing 
"conflicting results" from earlier simulations in the Intro- 
duction the Authors significantly distort the literature by 
citing these earlier simulations as "strongly imbalanced" . 
The Authors themselves presented simulations with en- 
ergy imbalance of around 3, while 1^ presented data for 
energy imbalance down to 1.35. 

Alignment controversy - The Authors define their 
alignment by Eq. 4. Why this specific function has to 
be relevant for the spectral slope is unclear. The Au- 
thors previously claimed in [l5| that this alignment mea- 
sure reduces the energy transfer function {6w^^{6w^)^) 
from [16], but this was just a plausibility argument equiv- 
alent to saying that alignment essentially causes anti- 
correlation of reducing the above function. Similar 
plausibility arguments can be applied to other alignment 
measures, such as IM = (ISjwt)'^ — S{w^)'^\) / {5{w^)'^ + 
reported in [s, 17|. However, IM does not show 
such a high scale-depencency as the Author's alignment. 
Why their alignment measure is so special, compared 
to other measures, apart from showing approximate Z^/^ 
scaling is never discussed by the Authors. 
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